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' Abstract: We investigate properties of the most general "PT-symmetric non-Hermitian 
Ph' 

Hamiltonian of cubic order in the annihilation and creation operators as a ten parame- 

' ter family. For various choices of the parameters we systematically construct an exact 

^ . expression for a metric operator and an isospectral Hermitian counterpart in the same 

O^. similarity class by exploiting the isomorphism between operator and Moyal products. We 
elaborate on the subtleties of this approach. For special choices of the ten parameters 

^ ^ the Hamiltonian reduces to various models previously studied, such as to the complex cu- 

' bic potential, the so-called Swanson Hamiltonian or the transformed version of the from 

' below unbounded quartic — x^-potential. In addition, it also reduces to various models 

■ not considered in the present context, namely the single site lattice Reggeon model and 

QQ I a transformed version of the massive sextic ±x^-potential, which plays an important role 

. as a toy model to identify theories with vanishing cosmological constant. 

o 

> 

^ ' 1. Introduction 



Non-Hermitian Hamiltonians are usually interpreted as effective Hamiltonians associated 
with dissipative systems when they possess a complex eigenvalue spectrum. However, from 
time to time also non-Hermitian Hamiltonians whose spectra were believed to be real have 
emerged sporadically in the literature, e.g. the lattice version of Reggeon field theory |||, |2[. 
Restricting this model to a single site leads to a potential very similar to the complex cubic 
potential V = ix^. Somewhat later it was found ||3| for the latter model that it possess 
a real spectrum on the real line. More recently the surprising discovery was made Q 
that in fact the entire infinite family of non-Hermitian Hamiltonians involving the complex 
potentials = z^{iz)'^ for n > possess a real spectrum, when its domain is appropriately 
continued to the complex plane. 

Thereafter it was understood [Q, ^ that the reality of the spectra can be explained by an 
unbroken T^T-symmetry, that is invariance of the Hamiltonian and its eigenfunctions under 
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a simultaneous parity transformation V and time reversal T. In case only the Hamiltonian 
is PT-symmetric the eigenvalues occur in complex conjugate pairs. In fact, the VT- 
operator is a specific example of an anti-linear operator for which such spectral properties 
have been established in a generic manner a long time ago by Wigner Q. However, in 
practical terms one is usually not in a position to know all eigenfunctions for a given non- 
Hermitian Hamiltoinian and therefore one has to resort to other methods to establish the 
reality of the spectrum. Since Hermitian Hamiltonians are guaranteed to have real spectra, 
one obvious method is to search for Hermitian counterparts in the same similarity class as 
the non-Hermitian one. This means one seeks similarity transformations r] of the form 

h = r]Hr]^^ = h) = rj^^H^^r] 4^ H^^ = rj^Hrj^'^, for rj = r]^ . (1.1) 



Non-Hermitian Hamiltonians H respecting the property (1.1) are referred to as pseudo- 
Hermitian Besides these spectral properties it is also understood how to formulate a 
consistent quantum mechanical description for such non-Hermitian Hamiltonian systems 
||8|, ^, ^] by demanding the ??^-operator to be Hermitian and positive-definite, such that it 
can be interpreted as a metric to define the r/-inner product. A special case of this is the 
CPT-inner product |^], which results by taking r/^ = CV with C = ^ (</>„|. For some 
recent reviews on pseudo Hermitian Hamiltonians see |l^, 11, |l^, |l4|. 



Since the metric-operator t/^ is of central importance many attempts have been made 
to construct it when given only a non-Hermitian Hamiltonian. However, so far one has 
only succeeded to compute exact expressions for the metric and isospectral partners in very 
few cases. Of course when the entire spectrum is known this task is straightforward, even 
though one might not always succeed to carry out the sum over all eigenfunctions. However, 
this is a very special setting as even in the most simple cases one usually does not have all 
the eigenfunctions at ones disposal and one has to resort to more pragmatic techniques, 



such as for instance perturbation theory ||1^, 16, 17, Rather than solving equations 
for operators, the entire problem simplifies considerably if one converts it into differential 
equations using Moyal products |jl^, 20, ^] or other types of techniques |2^]. Here we 



wish to pursue the former method for the most generic "PT-symmetric non-Hermitian 
Hamiltonian of cubic order in the creation and annihilation operators. 

We refer models for which the metric can be constructed exactly as solvable pseudo- 
Hermitian (SPH) systems. 

Our manuscript is organised as follows: In section 2 we introduce the model we wish to 
investigate in this manuscript, formulating it in terms of creation and annihilation operators 
and equivalently in terms of space and momentum operators. We comment on the reduction 
of the model to models previously studied. In section 3 we discuss in detail the method 
we are going to employ to solve the equations (p..l|), namely to exploit the isomorphism 
between products of operator valued functions and Moyal products of scalar functions. In 
section 4 we construct systematically various exact solutions for the metric operator and 
the Hermitian counterpart to H. As special cases of these general considerations we focus 
in section 5 and 6 on the single site lattice Reggeon model and the massive itx^-potential. 
In section 7 we provide a simple proof of the reality for the ix^"'"'"^-potentials and some of 
its generalizations. We state our conclusions in section 8. 
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2. A master Hamiltonian of cubic order 



The subject of our investigation is the most general T^T-symmetric Hamiltonian, which is 
maximally cubic in creation and annihilation operators ,a, respectively, 



He = Xio^ a+\2a^ +\zaa-\-\i-\-i{X^a^ -\-\Qa-\-\Ta^ +\%a^ a) a+\()a) aa+\iQaaa). (2.1) 

The Hamiltonian is a ten-parameter family with Aj G M. It is clear that this Hamiltonian 
is "PT-symmetric by employing the usual identification a = {ujx + ip)/\/2u; and a} = 
{ujx — ip) I \/2cJ with the operators in x-space x and p = —idx- The effect of a simultaneous 
parity transformation 7^ : x — > —x and time reversal T : t — > — t, i — > — i on the creation 
and annihilation operators is VT : a ^ —a^ ^ —c^ . Without loss of generality we may 
set the parameter uj to one in the following as it is simply an overall energy scale. 

In terms of the operators x and p the separation into a Hermitian and non-Hermitian 
part is somewhat more transparent and we may introduce in addition a coupling constant 
fl' G M in order to be able to treat the imaginary part as perturbation of a Hermitian 
operator. In terms of these operators the most general expression is, as to be expected, yet 
again a ten-parameter family 



t.t.t. 



{p ,x] {v,x] 
OL-j — hag — n h CKgx + Qiox 



(2.2) 
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Table 1: Special reductions of the Hamiltonian He- The map z{x) is defined in equation (6^) and 
g,a,A,m S M are coupling constants of the models. We abbreviated ki — — 4(16Ai + A2) and 
K2 = -4(48Ai + A2). 

We have symmetrized in He terms which contain p and x by introducing anticommu- 
tators, i.e. {^4, B} = AB + BA. This allows us to separate off conveniently the real and 
imaginary parts of H^ by defining 



Hcix,p) = ho{x,p) + ighi{x,p), 



(2.3) 



with /ig = ho and h\ = hi. In addition, the symmetrized version (2^) will lead to very 



simple expressions when we convert products of operator valued functions into expressions 
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involving scalar functions multiplied via Moyal products. For our definition of the Moyal 
product it implies that the parameters a do not need to be re-defined. Depending on the 
context, one (^]^) or the other (|2.2| ) formulation is more advantageous. Whereas the usage 
of creation and annihilation operators is more proned for an algebraic generalization, see 
e.g. |2^, the formulation in terms of operators x and p is more suitable for a treatment 
with Moyal brackets. The relation between the two versions is easily computed from 
the aforementioned identifications between the a,a^ and x,j5 via the relations a= MX and 
A = M~^a with M being a 10 x 10-matrix. Below we will impose some constraints on the 
coefhcients a and it is therefore useful to have an explicit expression for M at our disposal 
in order to see how these constraints affect the expression for the Hamiltonian in ( p.ll) . We 
compute the matrix 
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(2.5) 



which also exists for g ^ since det M = —g~^. 

The Hamiltonian H(. encompasses many models and for specific choices of some of 
the Qj it reduces to various well studied examples, such as the simple massive ix-potential 
1 24] or its massless version, the so-called Swanson Hamiltonian [p5| , p^ , 27, 21, the 
complex cubic potential together with his massive version and also the transformed 
version of the — x^-potential [^]. As we will show below, in addition it includes several 
interesting new models, such as the single site lattice version of Reggeon field theory p9|. 
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which is a thirty year old model but has not been considered in the current context and the 
transformed version of the ibx^- potential, which serves as a toy model to identify theories 
with vanishing cosmological constant [^0[. The latter models have not been solved so far 
with regard to their metric operators and isospectral partners. Besides these models, 
also includes many new models not considered so far, some of which are even SPH. 

To enable easy reference we summarize the various choices in table 1. 

Most SPH-models which have been constructed so far are rather trivial, such as the 
massive ix-potential or the so-called Swanson Hamiltonian. The latter model can be ob- 
tained simply from the standard harmonic oscillator by means of a Bogolyubov transfor- 
mation and a subsequent similarity transformation, which is bilinear in a and . Beyond 
these maximally quadratic models, the complex cubic potential was the first model which 
has been studied in more detail. Unfortunately so far it can only be treated perturbatively. 
The transformed version of the from below unbounded -z^-potential is the first SPH-model 
containing at least one cubic term. Here we enlarge this class of models. As a special case 



we shall also investigate the single site lattice version of Reggeon field theory |29| in more 
detail. Before treating these specific models let us investigate first the Hamiltonian in 
a very generic manner. 

Our objective is to solve equation ( |1.1| ) and find an exact expression for the positive- 
definite metric operator rf., subsequently to solve for the similarity transformation and 
construct Hermitian isospectral partner Hamiltonians. 

3. Pseudo-Hermitian Hamiltonians from Moyal products 
3.1 Generalities 

Taking solely a non-Hermitian Hamiltonian as a starting point, there is of course not a one- 
to-one correspondence to one specific Hermitian Hamiltonian counterpart. The conjugation 
relation in (^]^) admits obviously a whole family of solutions. In order to construct these 
solutions we will not use commutation relations involving operators, but instead we will 
exploit the isomorphism between operator valued function in x and p and scalar functions 
multiplied by Moyal products in monomial of scalars x and p. We associate to two arbitrary 
operator valued functions F{x,p) and G{x,p) two scalar functions F{x,p), G{x,p) £ S such 
that 

F{x,p)Gix,p) ^ F{x,p)i.G{x,p), (3.1) 
where S is the space of complex valued integrable functions. Here we use the following 



standard definition of the Moyal product ★ : 5 x 5 — > 5, see e.g. [31, 32, 21 1, 

F{x,p) * G{x,p) = F(3;,p)e5('^-^f-^p^-)G(x,p) (3.2) 
= f;^^^E(-l)* ( J did;''F{x,p)d^-'dlG{x,p). 

The Moyal product is a distributive and associative map obeying the same Hermiticity 
properties as the operator valued functions on the right hand side of ( |3.1| ), that is {F -k 
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G)* = G* -k F* . Following standard arguments we provide now explicit representations for 
the F(x,p) and F{x,p). We may formally Fourier expand an arbitrary operator valued 
functions F{x,p) and scalar functions F{x,p) as 

/oo fOO 
dsdtf{s,t)e'^'^+'P^ and F{x,p) = dsdtf{s,t)e'^'''+'P\ (3.3) 
-oo J — oo 

respectively. In terms of this representation the multiplication of two operator valued 
functions yields 

/oo 
dsdtds'dt' f{s, t)fis', t')et(*^'-*''^)e'(^+^')*+^(*+*')P, (3.4) 
-oo 

which follows using the identities e*("^+*P) = e''^/^e'^Pe''^/^ and e''^/^e'*P = e'^Pe''^/^e''K 
Is now straightforward to verify that the definition of the Moyal product ( |3.2| ) guarantees 
that the isomorphism ( |3.lD holds, since F{x,p)-kG{x,p) yields formally the same expression 
as ( |3.4[) with x,p replaced by x,p. 

The Hermiticity property is important for our purposes. We find that 

F\x,p) = F{x,p) ^ F*{x,p) = F{x,p). (3.5) 



This is easily seen by computing F^{x,p) using the representation ( p.3p . Then this function 
is Hermitian if and only if the kernel satisfies /* {s, t) = /(— s, —t), which in turn implies that 
F{x,p) is real. Positive definiteness of an operator valued function F{x,p) is guaranteed if 
the logarithm of the operator is Hermitian, that is we need to ensure that log-F(x,p) is real. 
Furthermore, it is easy to see that F{x,p) is PT-symmetric if and only if f*{s, t) = f{s, —t). 

As an instructive example we consider F{x,p) = for which we compute the 

corresponding kernel as f{s,t) = i'^~^"^6^"^\s)6^"'\t). From this it is easy to see that 
(ij;)™p" is "PT-symmetric, since i"^f{s,t) satisfies [i'^f{s,t)]* = i"^f{s,—t). 

In the present context of studying non-Hermitian Hamiltonians this technique of ex- 
ploiting the isomorphism between Moyal products and operator products has been exploited 
by Scholtz and Geyer 20 1, who reproduced some previously known results and also in 



|21], where new solutions were constructed. In p9| , |20| a more asymmetrical definition than 
( |3^ ) of the Moyal product was employed, i.e. F{x,p) * G{x,p) = F{x,p)e^^'^^pG{x,p). 
In comparison with ( |3.2| ) this definition leads to some rather unappealing properties: i) 
the loss of the useful and natural Hermiticity relation, i.e. {F * G)* ^ G* * F* , ii) the 
right hand side of the isomorphism in (^]^) is replaced by the less transparent expression 
F*{x,p) = e~'^^'^^^F{xjp) and iii) in ||2^ it was shown that the definition * leads to more 
complicated differential equations than the definition -k. The representation for the opera- 
tor valued functions F{x,p), which satisfies the properties resulting from the definition * 
differs from ( |3.3D by replacing e^i'^^+^p) — > Qi-sx^ttp -^^ Fourier expansion. 

3.2 Construction of the metric operator and isospectral partners 



We briefly recapitulate the main steps of the procedure [19, |2^, 21] of how to find for a 
given non-Hermitian Hamiltonian H a metric operator r]'^{x,p), a similarity transformation 
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r]{x,p) and an Hermitian counterpart h{x,p) using Moyal products. First of all we need to 
solve the right hand side of the isomorphism 

H\x,p)rj^{x,p) = rj^{x,p)H{x,p) = H'^ {x^p)-krj^{x,p) = rj^{x,p)-kH{x,p) (3-6) 

for the "scalar metric function" •if'{x,p). Taking as a starting point the non-Hermitian 
Hamiltonian H{x,p), we have to transform this expression into a scalar function H{x,p) 
by replacing all occurring operator products with Moyal products. We can use this ex- 
pression to evaluate the right hand side of the isomorphism of ( |3.6| ), which is a differential 
equation for rf'{x,p) whose order is governed by the highest powers of x and p in H{x,p). 
Subsequently we may replace the function rf'{x,p) by the metric operator rf'{x,p) using 
the isomorphism ( p.!] ) now in reverse from the right to the left. Thereafter we solve the dif- 
ferential equation r]'^{x,p) = r]{x,p)-kr]{x,p) for r](x,p). Inverting this expression we obtain 
r]~^{x,p), such that we are equipped to compute directly the scalar function associated to 
the Hermitian counterpart by evaluating 

h{x,p) = rj{x,p) -k H{x,p) -k ri~^{x,p). (3.7) 

Finally we have to convert the function r]{x,p) into the operator valued function r](x,p) 
and the "Hermitian scalar function" h{x,p) into the Hamiltonian counterpart h{x,p), once 
more by solving ( ^.1| ) from the right to the left. 

So far we did not comment on whether the metric is a meaningful Hermitian and pos- 



itive operator. According to the isomorphism (3.5) we simply have to verify that rf'ix^p), 
r]{x,p) and h{x,p) are real functions in order to establish that the corresponding operator 
valued functions rf{x,p), ri{x,p) and h{x,p) are Hermitian. We may establish positive 
definiteness of these operators by verifying that their logarithms are real. 

3.3 Ambiguities in the solution 

Obviously when having a non-Hermitian Hamiltonian as the sole starting point there is 
not a unique Hermitian counterpart in the same similarity class associated to the adjoint 
action of one unique operator rj. Consequently also the metric operator r/^ is not unique. 
The latter was pointed out for instance in ||2^ and exemplified in detail for the concrete 
example of the so-called Swanson Hamiltonian in |2^. In fact, it is trivial to see that any 
two non-equivalent metric operators, say rj"^ and r)^, can be used to construct a non-unitary 
symmetry operator S := rj^'^'q'^ = ffr]'"^ for the non-Hermitian Hamiltonian H 

= ifHrj-^ = fHf]-^ ^ [S, H] = [S^H^ = 0. (3.8) 

We may solve ( |3.8D and express one metric in terms of the other as 

ff' = (h'^Y rfH"" fornGN. (3.9) 



Thus we encounter here an infinite amount of new solutions. Likewise this ambiguity can 
be related to the non-equivalent Hermitian counterparts 

h = rjHr}-^ ,h = fjHfj-^ ^ [s,h] = [s,h] = 0, (3.10) 
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with symmetry operators s = rjVj^'^ri and s = T)^^ri^r]^^ . When rj^ = r] and r)^ = r) we 
obviously also have = s and = s. The expression for the symmetry operator s for h 
was also identified in 1 32 ] . 

There are various ways to select a unique solution. One possibility Q is to specify one 
more observable in the non-Hermitian system. However, this argument is very impractical 
as one does not know a priori which variables constitute observables. 

4. SPH-models of cubic order 

Let us study Hc{x,p) by converting it first into a scalar function Hc{x,p). Most terms are 
non problematic and we can simply substitute x ^ x,p ^ p, but according to our definition 
of the Moyal bracket ( |3.2| ) we have to replace p'^x p^-kx = p^x—ip, px^ — > p-kx"^ = px'^—ix, 
px ^ p -k X = px — i/2 etc. Replacing all operator products in this way we convert the 
Hamiltonian Hc{x,p) in ( |2.2[ ) into the scalar function 

Hc{x,p) = aip^ + a2p'^ + a3px'^ + a4p + a5x'^ + aQ + ig{a7p'^x + aspx + agx'^ + aiox). (4.1) 

Substituting ( [4.1| ) into the right hand side of the isomorphism into (|3.6| ) yields the third 
order differential equation 



aspxdp + a^xdp + ^d^dp + - a2pdx - ^aip^d-^ - ^x^d^ - ) rj^ (4.2 



agx + aiQX + aspx + ajp x + —pdxdp + —dxdp — -^xd^. — -agxSpJ rj 

for the "metric scalar function" ri'^{x,p). There are various simplifications one can make 
at this stage. First of all we could assume that either x or p is an observable in the non- 
Hermitian system, such that r]'^{x,p) does not depend on p or x, respectively. As pointed 
out before it is not clear at this stage if any of these choices is consistent. However, any 
particular choice p or x will be vindicated if (|4.2| ) can be solved subsequently for ri'^{p) or 
rf'{x), respectively. Here we will assume that rf'{x,p) admits a perturbative expansion. 
Making a very generic exponential "PT-symmetric ansatz, which is real and cubic in its 
argument for rf'{x,p) = e'xp g{qip^ + q2px'^ + QzP^ + liX^ + q^p) , we construct systematically 
all exact solutions of this form. Substituting the ansatz into the differential equation ( [4.2| ) 
and reading off the coefficients in front of each monomial in x and p yields at each order 
in g ten equations, by solving these equations we find five qualitatively different types of 
exact solutions characterized by vanishing coefficients a, and some additional constraints. 
We will now present these solutions. 

4.1 Non- vanishing j5x^-term 
4.1.1 Constraints 1 

We consider the full Hamiltonian Hc{x,p) in ( |4.1| ) and impose as the only constraint that 
the px^-term does not vanish, i.e. as 7^ 0. For this situation we can solve the differential 
equation ([4.2|) exactly to all orders in perturbation theory for 



Hc{x,p) = ho[x,p) + ig[ p X -\ px + uqx -\ x), (4.3) 

as "3 03 
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where we imposed the additional constraints 

aiag = 0307, 0209 = Oi50t^ + 0303 and 0409 = a^a% + 03010- (4-4) 



In (4.3) we have replaced the constants 07,08 and oio using ( |4.4| ). The solution of the 



differential equation is the metric scalar function 

r?2(^^^) = e-KiP'+^^+t^''). (4.5) 



Since rf'{x,p) is real it follows from ( |3.5D that the corresponding metric operator is Hermi- 
tian. Next we solve r]{x,p) -kr]{x,p) = T]'^{x,p) for r]{x,p). Up to order g'^ we find 



, , , 07p^ + asP + x'^ag 9 / ag («7 + ^a-jp^x"^ + 2aspx'^) + Og 
r][x,p) = I- g hff^ 



203 ■ " V 



3 



2^2 

3 



The corresponding Hermitian counterpart corresponding to this solution is computed by 
means of (13.71) to 



, , s 2 , 2 , , "3a7 3 , («5"7 + a3"8) 2 , ("508 + a3aio) 

hc{x,p) = a3px +05X +06H p H p H p 

09 09 09 

2(2 07P +03) (p (07^ +08) + 09^2 +010) , 4s 

-5 z + 0[g). (4.7) 

403 

Notice that since we demanded 03 to be non-vanishing these solutions can not be reduced 
to any of the well studied models presented in table 1, but represent new types of solutions. 
We may simplify the above Hamiltonians by setting various os to zero. 

Demanding for instance that x is an observable in the non-Hermitian system we are 
forced by (4^) to set 07 = og = and by ( [4. 41 ) also 01 = 02 = 0. The Hamiltonian in ( ^I^ ) 



then simplifies to 

TT / \ 2 2 • / 3 O4O9 , / A n\ 

Hc[x,p) = a^px + 04P + 05X + Og + ^^(oga; H x). (4.8) 

03 

Since ri'^{x,p) only depends on x in this case, we can compute exactly ri{x,p) = e ^^03^ ^ 
The Hermitian counterpart results to 

hc{x,p) = hQ{x,p) = a^px^ + a^p + osx^ + og. (4-9) 

If we require on the other hand that p is an observable, we have to choose 09 0. 
However, in that case the constraints (^^) imply that the non-Hermitian part of the Hamil- 
tonian ( |4.3| ) vanishes, i.e. we obtain the trivial case Hc{x,p) = /io(x,p). 
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4.1.2 Constraints 2 

In the construction of the previous solution some coefficients had to satisfy a quadratic 
equations in the parameters to guarantee the vanishing of the perturbative expansion. 
The other solution for this equation leads to the constraints ai = ay = 0, such that the 
non-Hermitian Hamiltonian simplifies. If we now impose the additional constraints 

tt3«io = onag and asas = 2a2a9, (4-10) 
we can solve the differential equation (|4.2|) exactly. For 



TT / \ 2 1 2 ■ ,2a2a9 o a4a9 , , , 

Hc(x,p) = a2P + aspx +aAp + a^x +aQ + ig[ px + agx H x) (4.11J 

as as 

we compute the exact scalar metric function to 



r/2(x,p) = e ^"3'''. (4.12) 

Clearly ri'^(x,p) is a Hermitian and positive definite operator, which follows from the facts 
that rf'{x^p) and logr/^(x,p) are real, respectively. Notice the fact that the Hamiltonian 
(4.11) does not follow as a specialization of (4.3), since the constraints ( [4.10| ) do not result 
as a particular case of ([4.4|). The Hermitian Hamiltonian counterpart corresponding to 
( 4.11 ) is computed with rj^{x,p) = e ^^"a^ by means of ( ^?7| ) to 

2 

1 / \ 9 9 9 9 ^2 9 / \ 

hc(x,pj = a2P + a^px + a^p + a^x + + g ^x. (4.13) 

a3 

Once again we may simplify the above Hamiltonians by setting various as to zero or other 
special values, except for the case ag — > for which the constraints (4.10) reduce the non- 
Hermitian part of the Hamiltonian ( [4.11| ) to zero. 
Thus this case requires a separate consideration: 

4.2 Non-vanishing px^-term and vanishing x'^-term 

Let us therefore embark on the treatment of the complementary case to the previous 
subsection, namely as 7^ and ag = 0. For these constraints we can solve the differential 
equation (4.2) exactly for the Hamiltonian 

Hc{x,p) = ho{x,p) + ig{a^p'^x + a^px + "^^"^"^ "^"'''^ x), (4.14) 

"3 

when we impose one additional constraint 

aioos = "5(0308 - a^a-j). (4.15) 

The "metric scalar function" results to 

7?2(x,p) = r]'^{p) = e^V^^ + ^ P), (4.16) 

Once again rj^{x,p) is a Hermitian and positive definite operator, which follows again 
from the facts that rj'^{x,p) and logr/^(x,p) are real. Since rf{x,p) only depends on p, 
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we can simply take the square root to compute r]{p). Then the corresponding Hermitian 
counterpart is computed by means of ( |3.7D to 



^^ = ^^+3 [j^P + P + P+ j 

(4.17) 



In fact we can implement the constraint (4.15) directly in the solution. The function 



rf{p) = {pa^ + as) ^ e^V^^ + ^ (4.I8) 

solves ( [4.2p for the generic Hamiltonian (^]^) with the only constraint that 03 7^ and 
og = 0. In this case the corresponding Hermitian counterpart is computed to 

4 [pa^ + 05) 



Implementing the constraint ( 4.15 ), the Hamiltonian ( 4.19 ) reduces to the one in ( 4.17| ). 



Similarly as the model of the previous subsection, these solutions can not be reduced to any 
of the well studied models presented in table 1, since 03 is assumed to be non- vanishing. 

4.3 Vanishing px^-term and non-vanishing x^-term 

Next we consider the complementary case to the previous two section, that is we take 
03 = in ( |2.2| ). For this set up we can only find an exact solution when we demand in 
addition that 05 7^ and ag = 0. For the non-Hermitian Hamiltonian 

Hc{x,p) = aip^ + a2p'^ + ctip + a^x^ + 00 + igiaip^x + a^px + oiqx ), (4.20) 

we can solve the differential equation ([4.2| ) exactly by 

r/2(x,p) = (3^^'+^?'"+^^') . (4.21) 



Once again 7]'^{x,p) only depends on p and we can simply take the square root to compute 
ry(p). Using ( |3.7D the corresponding Hermitian Hamiltonian is subsequently computed to 

hc{x,p) = aip + a2P +a4p + a5X + + g —. (4.22) 

4q!5 

Obviously these solutions can be reduced to various cases presented in table 1, notably the 
transformed — z^-potential and the Swanson Hamiltonian. 

4.4 Vanishing px^-term and non-vanishing p-term or non-vanishing p^-term 

Finally we consider the complementary case of the previous section by taking 03 = and 
allowing 05 to acquire any value. To be able to find an exact solution we need to impose 
the additional constraints 

CKi = Q7 = ag = 0, and a^ag = 2a2ai0) (4.23) 
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i.e. we consider the non-Hermitian Hamiltonian 



Hc{x,p) = a2P^ + aip + Q5X^ + a^+ig {a^px + aiQX ) , 
for which we can solve equation ( |4.2[) by 



rl^{x,p) 
rf{x,p) 



-gau)/oL4x'^ 



for 04 / 0, 
for a2 / 0. 



(4.24) 



(4.25) 
(4.26) 



As rf'{x,p) only depends on x, we can take the square root to compute r/(x) and subse- 
quently evaluate the corresponding Hermitian counterpart using (^ 



hc{x,p) 



2 2 2'^2Ctin 2 

a2P + a4P + a^x + ae + g 5 — x for 04 7^ 0, 

2 



hc(x,p) = a2P + 0:4^ + a^x + ae + g - — x for 02 7^ 0. 

4a2 



(4.27) 
(4.28) 



The Hamiltonian in ( 4.24| ) can be reduced to the Swanson Hamiltonian. Notice that 
when we impose ai = 04 = Q7 = aio = for the Hamiltonian in (4.20) and 04 = aio = 
for the Hamiltonian in ( [4.24 ), they become both identical to the Swanson Hamiltonian. 
The corresponding solutions for the metric operators reduce to r)^(x,p) = e^°*/^"^P and 
rf{x,p) = g-s"8/2o2a;^^ respectively, which are the well known non-equivalent solutions for 



the Swanson Hamiltonian, see e.g. |28]. This means according to (|3.8| ) we can identify a 
symmetry operator for the Swanson Hamiltonian as 



S{x,p) 



ae -,9 ao 



(4.29) 



Notice that S{x,p)-kH{x,p) = H{x,p)-kS{x,p) is a more difficult equation to solve in this 
example than ( |3.6D , since S is not of a simple exponential form as r/^. In fact this is what 
we expect. For instance supposing that the symmetry is some group of Lie type, a typical 
group element, when Gau/3 decomposed, would be a product of three exponentials. 



model\const 


ai 


02 


"3 


04 


05 


ae 




as 


ag 


aio 




ai 


02 




04 




ae 


aiag 




ag 


0409-0508 

OM 




4.11) 





02 




04 




ae 





2^209 


ag 


0409 
om 




4.14) 


ai 


OL2 




04 




ae 


ay 


as 





05(0308-0507) 






a\ 


02 





Q!4 




ae 


ay 


as 





aio 


^( 


4.24) 





02 







as 


ae 





202010 
04 





aio 


^( 


4.24) 










04 


05 


ae 





as 





2a2 



Table 2: SPH-models of cubic order. 

Obviously using the relation A = M^^a we can also convert our solutions into expres- 
sions using creation and annihilation operators. 
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5. Lattice version of Reggeon field theory 



Having studied the Hamiltonian (2^), (|2.2| ) in a very generic manner let us return to our 
original motivation and focus on some special cases, which have hitherto not been dealt 
with in the literature. It has been argued for more than thirty years that the lattice versions 
of Reggeon field theory (LR) |l| 



(5.1) 



with al, being standard creation and annihilation operators and A,g,g G M, possess a 
real eigenvalue spectrum [?]. This assertion was made despite the fact that the Hamiltonian 
in (15. 1|) is non-Hermitian. 



5.1 The single site lattice Reggeon model 

Many features of lattice models can be understood with a finite, possibly small, number of 



sites. Thus reducing the Hamiltonian in (5J) to a single lattice site yields the harmonic 
oscillator perturbed by a cubic perturbation in a and a^, such that unlike for the so- 
called Swanson Hamiltonian its properties can not be understood in direct analogy to the 
harmonic oscillator. The resulting model, which we refer to as single site lattice Reggeon 
model (SSLR) reads 

Hsslr{o., a^) = Aa^a + iga^ -|- a^^ a. (5.2) 

Obviously -ffssLR is a special case of the generic cubic "PT-symmetric Hamiltonians (p.l|). 
One may find already in the old literature, e.g. |34], that the Hermitian conjugation of H 
as specified in (|1 . 1|) can be obtained by an adjoint action with the parity operator 



r = e'-^^''' = r,\ (5.3) 
This is easily seen by noting that V acts on the creation and annihilation operators as 

VaV = -a and Va^V = -a) . (5.4) 



However, since the corresponding ri = \JV is not a Hermitian operator, we can not construct 
a Hermitian counterpart as specified in ( |1 . 1]) from this solution for . Another operator 
bilinear in a and a^, which has the same effect as V in (5^) when acting adjointly on a,a^, 
is 



7? 



e 2 



(aa—a^a^ ) 



(5.5) 



Taking now the square root of (^^) yields a Hermitian operator t] and serves therefore 
potentially for our purposes. Noting that this f\ acts as 



r/ar) ^ = ia) and r/a^r) ' = «a, 
such that a corresponding Hermitian counterpart is trivially obtained as 



h = fiHssLRfl ^ = -Aaa' + 



ga [a + a'] . 



(5.6) 



(5.7) 
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Clearly the spectrum of this Hamiltonian is not bounded from below for A > 0, resulting 
essentially from the fact that ff' in not a legitimate metric operator as it is not positive 
definite. Nonetheless, in the context of Reggeon field theory there is a considerable interest 
in the regime A < 0, such that the above argument contributes to an old discussion. 

Ignoring whether the metric is positive-definite or Hermitain we have two alternative 
solutions to equation (1.1) involving -ffssLR and therefore we have simple examples for the 
symmetry operators S = Vff and s = VjVfj in ( |3.8| ) and ( p. 10 ), respectively. 

Let us now try to find a more meaningful metric operator by using Moyal brackets. To 
commence, we have to convert the version (5.2) of the Hamiltonian into one which depends 
on X and p instead of a and a) . Using the aforementioned relations yields 

-f^sSLR(x,p) = y(p2 + - 1) + i-^{x^ + p^x -2x + ip) . (5.8) 

When ignoring the last three terms, this Hamiltonian becomes the massive version of the 
complex cubic potential of Bender and Boettcher Q. To proceed with our analysis we have 
to change H[x,p) from a function depending on operators to a scalar function. In most 
terms in ( |5.8D we can simply replace operators by scalars using x ^ x, p ^ p, but care 
needs to be taken with the term p'^x p^ -k x = p^x — ip. The resulting scalar function is 

Hs^hK{x,p) = a) a + iga) ★ + a)^ -ka = -^{'^^ — 1) + ^9^^^ ^-p^x — 2x), (5.9) 

where we have also scaled g —>■ g\f2 and set A = 1. For these values of the constants the 
differential equation (^) reduces to 

(2x6»p - 2pd^) r\^{x,p) = g (4x^ - 8x + Vx + 2pd.j,dp - 3x0^ - xd^) r\^{x,p). (5.10) 

Due to its close resemblance to the cubic potential we do not expect this equation to be 
exactly solvable and therefore resort to perturbation theory. 

5.2 SPH-models related to the single site lattice Reggeon model 

However, there are various models closely related to -ffssLR, which fit into the scheme of the 
previous section and are SPH. For instance, we can identify the exact solution of section 
4.1.1. by matching Hc{x,p) in (|4.3|) with Hsslr 



Hrg{a,a\ A, g, X) = HssLR{a,a\ A/2, g) +HssLR{-a,a\ -A/2, -X). (5.11) 



We find that all the constraints for the ten parameters in ( |4.4| ) are satisfied for this 
combination. An exact solution for the scalar metric function can then be identified as 

7/2(x,p) = e-»/^(^'+^''). (5.12) 



The Hamiltonian in ( 5.11 ) exhibits an interesting strong- weak symmetry 

F(y)(a,at,A,g,A) = i7Q(-a,at,-A,-A,-5). (5.13) 
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A further example for a SPH-model related to -f/^ssLR is 

HK2^{x,p,^,g,X) = HssLR(.x,p,A,g) - igx^, (5.14) 



which is identical to Hc{x,p) in ( ^.20| ) for certain values of the ten parameters Oj. Reading 
off those parameters yields as exact solution 

r]\x,p)=e-^3/Hf/3-2f)^ (5.15) 



for the metric operator according to ( [4.21| ). In fact, this model can be matched with 
the transformed version of the — z^-potential, for which the exact metric operator was 
constructed in 



5.3 Perturbative solution 

In order to solve the differential equation ( ^.lOl) perturbatively we make now the ansatz 

oo 

7/2(x,p) = 2j^5"c„(x,p) (5.16) 

n=0 

and the equation ( |5.10| ) is changed into a recursive equation for the coefficients 

{2xdp — 2pdx) Cn{x,p) = (4x^ — 8x + Ap^x + 2pdxdp — 3xOp — xd^) c„_i(x,p). (5.17) 

We may solve this successively order by order. Using the fact that limg_>o ?7^(x,p) = the 
initial condition is taken to be co(x,p) = 1. We then obtain recursively order by order 

ci{x,p) = p"^ — 2p + px'^ , (5.18) 
C2{x,p) = p^ — + + — 4p^x^ + 2p^x^ + p'^x'^, 

C3{x,p) = -p^ - Ap'^ - 5p^ + 24p^ - 4p + Spx^ - Qp^x"^ - Sp^x"^ + 2p'^x'^ - px'^ + -p^x^ 
3 3 

-4/x*^ + 2p^x'^, 

1 8 

ca{x,p) = -/^ - - 12/ + 76/ - 5/ - 72/ + 24x^ - IS/x^ + 104/x^ 
3 3 

-28p^x^ - 8/x2 - 13x^ + 28p^x^ - 20p^x^ - 8/x^ + 2/x^ - Ap^x^ 

3 3 3 3 

Apart from the usual ambiguities, which were discussed in section 3.3, there are new 
uncertainties entering through this solution procedure. Apart from the different choice 
for the integration constants, it is evident that on the left hand side of ( 5.17] ) we can 
always add to c„(x,p) any function of the Hermitian part of HssLRix,p), i.e. Cn{x,p) — > 
Cn{x,p) + \{x^ + p^ — \) is also a solution of the left hand side of ( 5.17] ). We fix this 
ambiguity by demanding rf'{x,p,g) •krf{x,p, —g) = 1 for reasons explained in pT] ]. 

Next we solve the differential equation rj{x,p) -kr]{x,p) = 7]'^{x,p) for r]{x,p) by making 
the ansatz 

oo 

rj{x,p) = l + Y,9''qn{x,p). (5.19) 

n=l 



-15- 



...generic cubic VT -symmetric non-Hermitian Hamiltonian 



We find order by order 

qi{x,p) = ci{x,p), q2{x,p) = C2{x,p)/2, (5.20) 

q3{x,p) = -p-' -p' - — + IQp" -3p- —p-^x +-p'x + —— + — -px 

D 4 i z 2 o 4 

-p^x'^ + Upx'^ - 2p^x^ 

, , 2 11 4 51 6 9 8 1 10 1 12 25p'^X^ 39 2 4 1 4 8 
qi{x,p) = -35/ + + —p^ - -p^ - + — ^ _^2^4 ^ _p4^8 

+ lpl0^2 _ 61 4 _ 23 4 4 _ 25 6^2 _ 7 2 6 _ 1 4^6 + 1 6^6 

6 8 2 ^ 2^ 4 2^ 3 6 

+45/x^ - p^x^ - p^x^ + 13x^. 

We are now in the position to compute the Hermitian counterpart to HsshR{x,p) by means 

of Q 

hssLRix,p) = ^ix^+p^- l)+g^ 0/ - 4p2 + 1 _ 4x2 ^ 3^2^2 ^ 3^4^ (^ ^l) 

("IZ/ _ 34p4 + 4p2 + 8 + 4x2 - 48p2x2 + lip^^^ _ ^4^4 ^ 31^2^4 ^ T^eA ^ ^(^6) 

Finally we recast our solution again in terms of creation and annihilation operators. Up to 
order the square root of the metric becomes 

7] = I + iV2ga^{a^ - a)a + g'^a^ a\2a^ a - - aa + 5)a - 2a) a) - 2aa + 2 a (5.22) 

and the Hermitian counterpart to the non-Hermitian Hamiltonian -ffsSLR acquires the form 

^SSLR = a^a + g^a^Ga^a + 4)a + / a^a^lOa^a'' + Waa - A8a^a)aa (5.23) 

+ a\20a^a^ + 20aa - l20aU)a - 32a^a + 0{f). 

As for all previously constructed perturbative solutions, it would be highly desirable to 
investigate in more detail the convergence properties of these solutions. 

6. Potentials leading to zero cosmological constants and the SSLR-model 

For the purpose of identifying vacuum solutions with zero cosmological constant 't Hooft 
and Nobbenhuis proposed in [^0| an interesting complex space-time symmetry transforma- 
tion between de-Sitter and anti-de-Sitter space 

dS adS : x^ ix^ = x ^ ix,p ^ —ip . (6.1) 

Since this transformation relates vacuum solutions with positive cosmological constant to 
those with negative cosmological constant, it can only be a symmetry for the vacuum if 
the cosmological constant is vanishing. In order to match this with a quantum mechanical 
Hamiltonian one demands that the map dS^adS sends H to —H, such that the vacuum 
state is the only invariant state of theory. This means any Hamiltonian of the form 

»rr, „ / - . ( rij + rrij = Akj form,- odd 

HAs{x,p) = y^,x''^f^^fj{x,p) with \ ' J r ' , 6.2 

^-^3 l^i + "^j = ^Kj for nij even, kj odd. 
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where the fj{x,p) = fj{ix,—ip) are arbitrary functions, is respecting this symmetry. A 
simple example for such a Hamiltonian was proposed in this context by Jackiw, see [3C], 



Hj{z, n, Ai, A2) = ^pI + Aii^ + X2z\ 



(6.3) 



Clearly dS— >adS maps Hj —Hj. In fact, his Hamiltonian involving a sextic potential 
was investigated before in 36, Notice that Hj is also 'PT-symmetric. Further- 
more, as pointed out in ^ for 1, 2, — 3) the groundstate wavefunction acquires 
a very simple form ipQ = exp(— z''/2) and Hj factorizes, such that it can be interpreted 
as the bosonic part of a supersymmetric pair of Hamiltonians. Moreover, this model is 
quasi-exactly solvable, meaning that a finite portion of the corresponding eigensystem has 
been constructed Bender:1996at. As discussed in Q one can continue the Schrodinger 
equation away from the real axis. Assuming an exponential fall off at infinity for Hj one 
may choose any parameterization which remains asymptotically inside the two wedges 



7 a 5 
--vr < 6 < vr 



and 



3 a 1 



(6.4) 



In fact, we can employ the same transformation as the one which was used successful for 
the — z'^-potential in |27] 

z{x) = -2iVl + ix. (6.5) 



For large positive x we find z ~ e ^'^/^ G Wj? and likewise for large negative x we find 
We then find 

Hi(x,p^,g) = Hi[z{x),p^,g\ = — ^^K^) = Y + T 2^ - 2 + 2^^^' ~ ^^'^^ 

This allows us to interpret the Hamiltonian Hj{z, 0,, Ai, A2) as a perturbation of the exactly 
solvable model H4{z,pz, 9), since 



Hj [z{x), 1,5/384,5/8] = Hi{x,p^,g) + f - |. 

6 



(6.7) 



We can also relate the special case to an exactly solvable model 



fx Px 



Hj[zix),l,2,-3] = ^-^ + 384x2- 116 + i ^ 
= Hf^{x,p)+il28x^ 



{x,pI} 



+ 372x + 128x'^ 



(6.8) 
(6.9) 



It would be very interesting to investigate also the possibility to have a p-dependence 
in the potential Hds{x,p), but unfortunately this always leads to equations with order 
greater than three and is therefore beyond our generic treatment. 
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7. A simple reality proof for the spectrum of + z'^{iz 



,2m+l 



Considerable efforts have been made to prove the reahty of the spectrum for the family of 
Hamiltonians Hn = + z'^{iz)^ for n > 0, see for instance ^]. Unfortunately most 
of the proofs are rather cumbersome and not particularly transparent. The simplest way 
to establish the reality of the spectrum for a non-Hermitain Hamiltonian is to construct 
exactly the similarity transformation, which relates it to a Hermitian counterpart in the 
same similarity class. So far this could only be achieved for the case n = 2 in a remarkably 
simple manner ||27|] . Here we present a trivial argument, which establishes the reality for a 
subclass of the massive version of Hn with n being odd. Considering 

H,^ = ^{p^+x^)+gx\ixf"'-\ A,5GM,meN, (7.1) 

2m+l/2 V / 

it is trivial to see that the metric operator r] = eT('^'^~'^^'^^) introduced in (|5.5|) transforms 
Hm adjointly into a Hermitian Hamiltonian 

hm = f,H"^r' = -Aaat + (a + at)'"""' = hi (7.2) 

= -^(j52 + x2)-5x^-+i. 

We can immediately apply the same argumentation to generalizations of the single site 
lattice Reggeon Hamiltonian. The non-Hermitian Hamiltonian 

H^SLR = Aata + |at (a + at) a (7.3) 
is transformed to the Hermitian Hamiltonian 

^SSLR = vHEhLRT' = -Aaat - ^^at (a + at)'"'^' a. (7.4) 

For the reasons mentioned in section 5 the operator fi'^ is not a proper positive-definite 
metric, but for the purpose of establishing the reality of the spectrum that is not important. 

8. Conclusions 

We have systematically constructed all exact solutions for the metric operator which is 
of exponential form with PT-symmetric real and cubic argument and adjointly complex 



conjugates the most generic Hamiltonian of cubic order ( |2.1| ), (2.2). Our solutions are 
characterized by various constraints on the ten parameters in the model. Several of the 
SPH-models may be reduced to previously studied models, but some correspond to en- 
tirely new examples for SPH-models. We used the metric to construct the corresponding 
similarity transformation and its Hermitian counterparts. 

We have demonstrated that exploiting the isomorphism between operator and Moyal 
products allows to convert the operator identities into manageable differential equations. 
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Even when no obvious exact solution exists, perturbation theory can be carried out on the 
level of the differential equation to almost any desired order. An important open question, 
which will be left for future investigation is concerning the convergence of the perturbative 
series. There are obvious limitations for the method as it works only well for potentials of 
polynomial form, as otherwise the differential equations will be of infinite order. In addition 
even for the differential equations of finite order not all solutions of all possible have been 
obtained as one has to make various assumptions. For instance the solution ( 4.1^ ) would be 



missed for an ansatz for the metric of a purely exponential form. Therfore it would be very 
interesting to compare the method used in this manuscript with alternative techniques, 



such as the one proposed in p^ . 



The special cases considered, namely the single site lattice Reggeon model (|5.2D as 
well as the sextic potential ( |6.3| ) are both not SPH within our framework. However, both 
of them may be understood as quasi exactly solvable models perturbed by some complex 
cubic potential. 

There are some obvious generalizations of the presented analysis, such as for instance 
the study of generic quartic, quintic etc Hamiltonians. 
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